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Abstract 


A  new  model  for  laser-induced  desorption  is  proposed  based  on 
a  master  equation  for  the  photon  population  in  a  closely-coupled  sub¬ 
system  involving  the  relevant  degrees  of  freedom.  The  validity  of 
approximations  such  as  the  Born  and  Markovian  approximations  is 
discussed  in  terms  of  the  laser  pulse  duration  and  multiple  time 
scales  of  the  adspecies-surface  system.  The  desorption  rate  is 
numerically  calculated  from  the  photon  population  and  the  threshold 
number  of  photons  absorbed  by  the  adspecies.  The  effects  of  anharmonicity , 


dephasing  and  laser  detuning  on  the  desorption  rate  are  examined.  The 


mechanism  of  dephasing  and  the  possible  direct  and  indirect  desorption 


channels  are  discussed.. 
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I.  Introduction 


Laser-stimulated  surface  processes  (LSSP)  have  been  investigated  during 

the  past  several  years ,  due  both  to  their  academic  interest  and  industrial 

1  2  3 

potential.  Recent  progress  in  experimental  and  theoretical  *  studies  and 
4 

in  applications  of  LSSP  has  been  reported.  Laser  excitation  and/or  desorp¬ 
tion  of  adspecies  have  been  investigated  theoretically  by  a  variety  of  tech- 
niques,  including  harmonic,  anharmonic  and  Morse'’  potential  models  and  master 
equation  approaches, ^  where  attention  has  been  focused  on  the  excitation 
of  the  active  mode  and  the  population  of  its  vibrational  states.  The  desorp¬ 
tion  rate  is  usually  overestimated  by  the  harmonic  model^  whereas  it  is  under- 

7  8 

estimated  by  the  one-dimensional  Morse  potential  model.  ’ 

In  the  present  paper,  we  propose  a  new  model  in  which  all  the  degrees 
of  freedom  of  a  closely-coupled  subsystem  (including  the  active  mode  of  the 
adspecies)  are  treated  on  an  equal  footing.  The  desorption  rate  is  calculated 
from  the  photon  population  and  the  threshold  number  of  photons  absorbed  by  the 
adspecies  as  a  whole,  rather  than  from  the  average  excitation  of  the  active 
mode.  In  Section  II,  the  generalized  master  equation  is  developed  within 
the  Born  approximation,  and  a  reduced  master  (rate)  equation  including  the 
anharmonicity  of  the  adspecies  surface  potential  is  investigated  within  the 
Markovian  approximation.  Numerical  results  are  shown  in  Section  III,  where 
dephasing-enhanced  desorption  is  suggested.  In  Section  IV,  the  mechanism  of 
dephaslng  and  a  variety  of  desorption  channels  are  discussed. 

II.  Master  Equation 

Let  us  first  define  our  model  system,  which  is  particularly  appropriate 
for  admolecules  with  very  fast  intramolecular  relaxation.  In  the  frequency 
domain,  the  total  adspecles-surface  system  is  divided  into  two  subsystems, 
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and  S2»  where  consists  of  the  active  mode  (A)  and  the  strongly-coupled 
modes  (B)  within  the  adspecies,  and  S2  consists  of  the  remaining  lower  fre¬ 
quency  modes  serving  as  a  heat  bath  (see  Fig.  1).  From  the  concept  of  the 
energy-gap  law*  we  expect  Intramode  relaxation  within  to  be  much  faster 
than  intersystem  relaxation,  i.e.,  w^ere  relaxa¬ 

tion  time  between  the  1  and  j  modes.  When  the  system  is  irradiated  by  a  laser 
field,  will  be  excited  through  the  active  mode  and  its  coupling  to  the  B 
modes.  The  excitation  rate  of  S^,  after  a  short  time  of  Y^,  is  proportional 
to  the  absorption  cross  section  of  the  active  mode.  We  can  therefore  express 
the  Interaction  between  the  laser  field  and  as 

HSF(t)  "  v<tMa++a>*  (1) 

where  V(t)<r  yABu,E(t)cos(ut)  is  proportional  to  the  derivative  of  the  active¬ 
mode  dipole  moment,  y * ,  the  coupling  factor  between  A  and  B,  y  ,  and  the  laser 
field  amplitude,  E(t),  with  a  frequency  u>  which  is  near-resonant  to  that  of  the 
active  mode,  a^  and  a  are  the  ladder  operators  for  defined  in  photon 
space,  i.e.,  a+|n>  - Xv+I  |  n+l> ,  where  |n> is  a  basis  function  for  as  a  whole, 
rather  than  for  just  the  active  mode  which  reaches  steady  state  in  the  time 


The  vibrational  Hamiltonian  of  the  system  is 


H  -  Hg  +  Hg  +  Hg  +  HgF(t),  (2) 

where  Hg  couples  and  S2  (it  is  generally  time-independent).  The  density 

lg 

matrix  for  the  total  system  obeys  the  Liouville  equation 


,+ls+lsf 


)  p. 


dp 

dt 


(3) 
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Within  the  Born  approximation,  L_  +  Lg  +  LS  +  LS  ’  and  for  Markovian 

31  a2  1  2 

processes  where  the  characteristic  time  of  the  heat  bath  (S2>  is  much  shorter 
than  other  time  scales  of  the  system,  the  S2  variables  may  be  eliminated.  We 
can  then  work  with  a  reduced  density  matrix,  p  =TrB(p),  which  obeys  the  equation 

f  *  -i,LSl+VlVLSFre.  <*> 

where  the  effects  of  the  interaction  between  and  S2»  Lg  ,  are  now  reduced 

to  a  constant  frequency  shift,  Lg,  and  a  damping  factor,  L^.  We  now  introduce  the 

Markovian  approximation  for  the  relaxation  associated  with  Lg  (assumed  to  be 

time  independent),  but  retain  memory  effects  due  to  the  off-diagonal  matrix 

elements  associated  with  the  laser  excitation.  For  a  Hamiltonian  such  as 

given  by  Eq.(l),  we  obtain  a  generalized  master  equation  for  the  diagonal 

o 

matrix  element,  P  s  [Tr_p]  ,  as  (detailed  discussions  of  the  generalized 
n  d  nn 

16—18 

master  equation  are  available  in  the  literature  ) 


dP 


~dt  *  ’'l""+1>Pn-"Pnl 


+  1  dt'IH  , (t')-P „(t’)J  +  W  , (c’)-P  <t') ]) 

q  n,n— i  n-i  n  n,n+l  n+i  n 


(5) 


with  the  time-dependent  photon  absorption/emission  rates  given  by 


Wn,n+1(T)  “  2(n+1)  n(t)  SHt’)e~irTcos(an+1T)  (6) 

rT 

W  i  (t  )  “  2n  0(012(1*)  e  11  cos(A  t),  (7) 

n»n— i  n 

Here,  fl(t)  ■  V(t)  /cos(ut) ;  x  •  t-t*  ;i  *  w  -  u  *  A  -  2cm,  where  u  is 

in  inn 

the  transition  frequency  between  levels  m  and  n,  A  -  w^q-  4)  is  the  laser 
detuning  with  respect  to  the  fundamental  frequency  and  e  is  the  anharmonicity 
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and  r  *  y^/2  +  is  the  total  damping  associated  with  the  off-diagonal 

matrix  elements,  where  and  Y  denote  the  energy  and  (pure)  dephasing 

factors,  respectively.  We  have  used  the  rotating-wave  approximation  and 

the  dipole  transition  for  the  n  dependence  of  W  (such  that  m  =  n  t  1)  for 

m,n 

absorption  (m  <n)  and  emission  (m  >n).  Note  that  is  the  probability  of 

the  adspecies  as  a  whole  (S^  subsystem)  absorbing  n  photons,  in  which  the 

active-mode  state  is  given  by  its  steady-state  excitation,  i.e.,  n^(s.s.)  « 

YAB<’I/^u),  where  a  is  the  steady-state  absorption  cross  section. 

For  a  very  short  laser  pulse  with  duration  t  <  r-1,  the  transient  solu- 

P  * 

tion  of  Eq. (5)  contains  the  memory  effects  of  the  population  function,  and  the 

Markovian  approximation  may  not  be  applied  to  this  time-dependent  excitation. 

We  shall,  instead,  focus  on  the  case  of  a  long  pulse  excitation,  e.g.,  a  CO 

laser  pulse  with  tp  =  10  ns  which  is  much  longer  than  the  dephasing  time  r_1, 

such  that  the  population  function  P^(t)  and  Q(t)  are  slowly  varying  and  may  be 

factored  out  of  the  integrals  in  Eq.(5).  Employing  this  adiabatic,  Markovian 

approximation,  the  transition  rates  W^  ^  have  no  memory  effects,  although  they 

are  still  time  dependent  due  to  the  laser  temporal  profile  E(t),  and  a  simpler 

9  16 

master  equation  is  then  obtained:  * 


dP  n(P  -P  ) 

sr  -  y  <"+1>pn+rr,p„i +  2rn2<t)  ( — 1"1  11 


<n+1)(Pn+1-Pn) 


[ A-2e (n-1) ] 2+  r  2  (a-2en)2+  r  2 

Several  features  of  this  equation  are:  (i)  the  absorption  cross  section  for 
the  adspecies  (S^  subsystem)  is  saturated  at  a  higher  photon  population  due  to 
the  anharmonicity;  (ii)  the  first  two  terms  describe  the  actual  energy  flow  from 
to  the  surface  ( S2 ) 5  (iii)  the  pure  dephasing  factor  O^)  changes  only  the 
intramolecular  phase  without  changing  the  energy  populations;  (iv)  the  multi¬ 
phonon  relaxation  factor  O^)  is  generally  strongly  temperature  dependent  and 


(8) 
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r 


is  a  strongly  decreasing  function  of  the  energy  gap  between  and  S^; 

(v)  the  remaining  terms  containing  the  factor  r  indicate  that  the  effect 
of  the  applied  field  is  always  intimately  related  to  the  phase  dissipation 
or  broadening  mechanisms. 

Equation  (8)  may  be  referred  to  as  the  usual  rate  equation  in  which 

no  memory  effects  are  preserved  either  in  the  relaxation  or  the  excitation 

process.  For  weaker  dephaslng  systems  subject  to  a  short  pulse  excitation, 

the  "most  general"  form  of  the  generalized  master  equation  is  given  by  Eq.  (5) 

18 

in  which  the  relaxation  terms  are  replaced  by 


2  t  -(iu>  +r  )(t-t’) 

72  Re  I  f  df  <  Hl(t)Hl(t')>  -P  (t')}e  mn  mn 


where  *  (Ym+Yn)/2  +  ^2  t*ie  tota^-  surface-induced  damping.  Furthermore, 
for  a  short-pulsed  laser  with  high  intensity,  the  Born  approximation  may  not 
be  valid.  Improvements  can  be  made  by  performing  a  trace  over  the  photon(laser)- 
dressed  states  of  the  heat  bath  (S2) •  For  the  case  of  strong  coupling  between 
and  S^,  may  not  be  used  as  a  perturbation,  and  again  the  Born  approxi¬ 
mation  is  questionable.  To  go  beyond  this,  we  may  carry  out  a  unitary  trans¬ 
formation  on  with  subsequent  laser  interaction  with  the  phonon (surface) - 
2 

dressed  states.  These  improvements  allow  for  both  multiphonon  (relaxation) 
and  multiphoton  (excitation)  processes.  The  absorption  cross  section  for  the 
subsystem,  for  example,  will  be  related  to  the  laser  intensity  by  a  power 
law  o  «  1°  for  n-photon  processes.  We  also  note  that  a  nonlinear  power  law, 
o  *  Im  with  m<  1,  is  also  possible  by  anharmonic  saturation.^ 


Desorption  Rate 


For  a  first-order  desorption  process,  the  rate  constant  is  related  to 
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the  Inverse  of  the  mean  first-passage  time,  t,  and  is  given  by 


where  n*  is  the  desorption  threshold  number  of  photons  absorbed  by  the  adspecie 
Numerical  results  are  shown  in  Figs.  2  and  3.  From  Fig.  2  we  see  that  the 
desorption  rate  is  a  strongly-decreasing  function  of  the  anharmonicity . 
Dephasing-enhanced  desorption  is  seen  in  Fig.  3,  where  an  optimal  value  at 
:  1.8  is  found  for  a  maximal  desorption  rate.  This  is  realized  by  the 
fact  that  the  dephasing  tends  to  compensate  for  anharmonicity  "bottleneck" 
effects.  An  enhancement  factor  of  about  3  is  possible  if  we  tune  the  dephasing 
to  the  optimal  value  of  1-8,  in  comparison  with  -  4.  We  note  that  the 
dephasing-enhanced  effects  caused  by  the  nonlinear  behavior  are  absent  in 
a  harmonic  model,  where  is  always  a  decreasing  function  of  y 

In  Fig.  4  we  show  the  time  evolution  of  the  desorption  probability, 
defined  by 

PD(t)  "  ^  Pn(t)'  (11) 

n-n* 

Again,  we  see  the  nonlinear  behavior  of  the  dephasing  effects.  For  a  fixed 
time,  say  t  *  100  units,  PD  increases  when  increases  from  1  to  5,  but  it 
decreases  at  higher  values  of  Y 
IV.  Discussion 
A.  Mechanism  of  Dephasing 

As  indicated  in  the  previous  section,  the  dephasing  factor  changes 
only  the  phase  of  the  excited  subsystem  without  changing  its  energy  popula¬ 
tions.  In  the  language  of  gas-phase  scattering  theory,  this  can  be  referred 
to  as  an  "elastic"  scattering  process.  In  the  presence  of  a  laser  field, 
dephasing  can  be  viewed  as  energy  transfer  from  the  excited  subsystem  to  a 
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"virtual"  state  of  the  field  and  back  to  the  original  excited  state.  Within 
the  adspecies/surf ace  system  itself,  the  dephasing  is  due  to  one  or  more  of 
the  following:  *  ’  (i)  dephasing  of  the  active  dipole  of  the  adspecies 

(or  dephasing  of  the  laser  field  itself);  (ii)  fluctuation  of  the  conformation 
of  the  adspecies  due  to  the  thermal  energy  provided  by  either  the  initial  sur¬ 
face  temperature  or  laser-heating  effects;  (iii)  librational  and  rotational 
relaxation  within  the  vibrational  manifold  of  the  excited  subsystem;  (iv)  intra¬ 
molecular  mode-mode  coupling  within  and  intermolecular  coupling  between 
and  (v)  lateral  motion  or  migration-induced  elastic  collisions  between 
adspecies  and  with  the  substrate  surface;  (vi)  surface  phonon-dispersion- 
induced  level  broadening  of  the  vibrational  manifold. 

We  note  that  the  dephasing  factor  introduced  in  Section  II,  is 
formally  derived  from  the  Liouville  equation  within  the  Markovian  approximation, 
where  the  laser  field  amplitude  is  assumed  to  be  generated  by  a  coherent  source. 
However,  the  field  Itself  may  often  be  only  partially  coherent  and  hence  will 
cause  additional  dephasing  of  the  excited  subsystem.  To  demonstrate  this 
mathematically,  we  start  with  the  reduced  density  matrix  equation  of  motion 

for  Pm„  =  m  *  n> 

mn  B  mn 


dP 

mn 

dt 


i  [w  +  uj(t)  ]P  , 
mn  mn 


(12) 


where  we  have  introduced  a  stochastic  frequency  shift  of  the  system,  u>(t), 
caused  by  the  partially  coherent  behavior  of  the  laser  field.  By  assuming  a 
"white  noise,"  «a(t)>  *  0,  and  a  Markovian  correlation,  <w  (t)u>(t')>  =  Y25(t~tf)* 
we  obtain  an  ensemble-averaged  (over  the  stochastic  variable)  equation 


dP 

mn 

<  — : > 


dt 


-  [i«  +  (r  +  0]<P  >, 

mn  mn  i.  mn 


where  ^  is  the  laser-induced  dephasing. 
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(13) 


» 


t 
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B.  Mechanism  of  Desorption 

The  desorption  mechanism, which  may  involve  several  channels,  depending 

on  the  nature  of  the  adspecies  (e.g.,  physisorbed  or  chemisorbed,  adatom  or 

admolecule,  and  the  potential  of  the  adbond),  and  the  associated  energy  levels  of  the 

system  are  shown  in  Fig.  5.  The  desorption  may  result  from  direct  laser 

excitation  or  indirect  energy  transfer  processes.  The  details  of  the  desorp- 

3  4 

tion  channels  from  Type  I  to  Type  VI  have  been  discussed  previously  ’  . 

As  a  concluding  remark,  we  note  that  a  complete  description  of  laser-induced 
desorption  should  at  least  include  the  competing  processes,  e.g.,  migration- 
induced  desorption  and  re-adsorption  after  the  mean  first-passage  time.  Fur¬ 
thermore,  for  systems  with  more  than  one  active  mode,  desorption  may  occur 
via  a  two-stage  process:  the  initial  stage  of  absorbing  n*  photons  followed 
by  thermal  phonon-assisted  and/or  V-V  energy  transfer  among  the  excited 
active  modes.  This  type  of  process  is  under  investigation  in  our  laboratories. 
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Figure  Captions 

1.  Schematic  diagram  for  the  frequency  spectrum  of  the  total  system  Sj  +  S2> 
where  -  A  +  B ,  with  the  corresponding  relaxation  factors  indicated. 

2.  Anharmonicity  (t)  dependence  of  the  desorption  rate  in  arbitrary  units. 

The  parameters  used  are  <V  VV  a)  "  (0.1, 2. 0,0. 1,0.0) ,  with  a  laser  pulse 
duration  of  tp  -  10  ns,  where  fl(t)  -  ^sin^t/tp) .  y2,  fi()  and  a  are 

in  the  unit  of  10^  s  and  n*  “3  has  been  used  in  Eq.(10). 

3.  Phase-relaxation  (>2)  dependence  of  K^.  The  parameters  used  are  (Y^,n^,e,A,t  ) 
"  (0.5,0.1,1,3,10).  The  envelope  function  of  the  laser  pulse  is  the  same 

as  in  Fig.  2. 

4.  Time  dependence  of  the  desorption  probability  PQ(t),  where  n*-3  has  been 
used  in  this  calculation.  A  cw  laser  field  has  been  assumed,  and  the 
other  parameters  are  (y^  nQ,  e,  a)  *  (0.01,  0.1,  1.5,  3.0). 

5.  Schematic  diagrams  of  adspecies-surf ace  systems  and  the  associated  energy 
levels,  where  A,  B  and  M  represent  the  adspecies  (adatom  or  admoleculc) , 

C  represents  the  substrate  (or  bath  modes),  and  the  laser  radiation  is 
indicated  by  the  wiggly  lines.  Several  types  of  desorption  channels  are 
illustrated:  (I)  direct  desorption  via  active-mode  excitation,  (II)  direct 
desorption  via  the  quasi-continuum,  (III)  indirect  desorption  via  tunneling, 
(IV)  indirect  desorption  via  substrate  heating,  (V)  phonon-assisted 
desorption  and  (VI)  indirect  desorption  via  dynamics. 
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